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We discuss a new method to extract neutrino signals in low energy experiments. In this scheme the 
symmetric nature of most backgrounds allows for direct cancellation from data. The application of 
this technique to the Palo Verde reactor neutrino oscillation experiment allowed us to reduce the 
measurement errors on the anti-neutrino flux from ~ 20% to ~ 10%. We expect this method to 
substantially improve the data quality in future low background experiments such as KamLAND 
and LENS. 



PACS 14.60.Lm, 14.60. Pq, 29.85.+C 

I. INTRODUCTION 

Backgrounds are a major concern in low-energy neu- 
trino experiments where signals have low rates and are 
easily mimicked by other phenomena. Several types 
of coincidence schemes, specific to particular neutrino- 
induced processes, have been proposed to improve the 
signal-to-noise ratio. One classic example is the use of 
the inverse-/? decay process 

v c + P — > c + + n (1) 

in liquid scintillator in the discovery of neutrinos and 
many subsequent experiments Positrons deposit 

their energy in the scintillator and annihilate, yielding 
two 511 keV 7's. Neutrons are captured after thermaliza- 
tion, producing 7's. The two parts of the event are sepa- 
rated in time by a delay ranging from tens to hundreds of 
microseconds depending upon the nucleus on which the 
neutron captures. The use of similar time correlations 
has been proposed for solar neutrino detection || . 

It is often the case that experiments are still back- 
ground limited even when such coincidence schemes are 
adopted, particularly when the signal rate is very low 
and can not be varied. Using data from the Palo Verde 
neutrino oscillation experiment Q, we have found that 
most backgrounds have a peculiar symmetry in the en- 
ergy depositions between the two parts of an event that 
is not present in the neutrino signal. Such symmetry al- 
lows one to eliminate most of the background by direct 
subtraction with the data itself. 

In this paper we discuss in detail the method using 
the data from the Palo Verde experiment as an example 
and its application to future experiments such as Kam- 
LAND H and LENS || , where signal rates are expected 
to be substantially lower. 

II. BACKGROUNDS TO REACTOR NEUTRINO 
EXPERIMENTS 

Low-energy electron anti-neutrinos from nuclear reac- 
tors are unique tools to study oscillations in the regime 



of large mixing angle and small mass differences. Such 
combination of parameters has recently received a good 
deal of attention as it is consistent with a number of ob- 
servations involving solar and atmospheric neutrinos ||. 
In recent times two experiments of this type have 
been set up to search for v c — 9 X oscillations compat- 
ible with the atmospheric neutrino anomaly. In these 
experiments, electron anti-neutrinos from reactors with 
energies less than 10 MeV, are detected by the reaction 
(1) in liquid scintillator. There are two types of back- 
grounds to this process: uncorrelated background from 
environmental radioactivity randomly occurring during 
the time coincidence window, and correlated background 
from cosmic-muon-induced fast neutrons. While the first 
one can be easily measured by varying the time corre- 
lation window, the second one is more difficult to be 
measured unambiguously. Neutrons discussed here are 
produced cither in the laboratory walls or inside the de- 
tector. Michel electrons from muon decays are not a 
background since their time correlation is short and their 
energy deposition too large. 

Fast neutrons can mimic the anti-neutrino signal in the 
following two ways: 

• one-neutron background: A proton recoil is pro- 
duced through a fast-neutron scattering mimick- 
ing the e + signature; the neutron is then thermal- 
ized and captured like in the case of anti-neutrino 
events. 

• two-neutron background: The fast neutron can pro- 
duce a secondary neutron through a spallation pro- 
cess on nuclei; both neutrons arc then captured 
simulating the two parts of an anti- neutrino event. 

Both backgrounds are very difficult to measure except 
in the case when the v e source (in this case nuclear reac- 
tors) can be turned off, hence eliminating the signal. This 
favorable circumstance was available only to the Chooz 
experiment [0. Generally, theoretical models describing 
neutron production are not considered accurate enough 
to provide a viable tool for background subtraction. 
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Fast neutrons are produced mainly in muon capture 
and muon spallation. While the first process is well un- 
derstood, the second is poorly known. Although the to- 
tal neutron yield from muon spallation has been, to some 
extent, experimentally measured theoretical mod- 

els ]l0|-|l2"[ are not consistent with each other and with 
data. In addition the few measurements of the neutron 
energy spectrum ]l3|] a re not well reproduced by theo- 
retical calculations |l4||l5| ]. The interpretation of experi- 
mental data is complicated by the fact that the neutron 
energy spectrum depends upon the muon spectrum that, 
in turn, is a function of the depth at which the measure- 
ment was carried-on. 

Since modern reactor neutrino oscillation experiments 
typically have long baselines and observe anti-neutrinos 
from more than one reactor, in most cases it is impossi- 
ble to completely turn off the signal source. Hence many 
experiments rely on the power variations that generally 
occur during refueling of some of the reactors in order 
to subtract backgrounds. This method, here referred to 
as "ON-OFF method" , has serious limitations since: a) 
the statistical error is large as only a small fraction of 
the neutrinos are used as signal while most of them are 
subtracted away with the background; the smaller the 
power excursion the larger the statistical error; b) since 
reactors are kept to full power for a very large fraction of 
time (because of obvious economic reasons) statistical er- 
rors are dominated by the short low power periods while 
the majority of the data taken by the experiment is not 
useful to improve the measurement accuracy; c) the sub- 
traction method only works under the assumption that 
backgrounds are stable over the periods of several months 
that separate the full-power periods from the low-power 
ones; d) complete systematic checks on data can only be 
done after an entire reactor cycle that generally corre- 
sponds to a period of six months to one year. 

The new technique, that we call the "swap method" , 
avoids such limitations. 

III. THE SWAP METHOD 

The swap method uses symmetries of the data to di- 
rectly eliminate most of the backgrounds and a Monte 
Carlo calculation to estimate the residual background. 
The same symmetries that guarantee the cancellation in 
data also make the whole process rather insensitive to 
imperfections of the Monte Carlo model. 

We first select neutrino events by requiring the prompt 
part as positron-like and the delayed part as neutron-like. 
We have: 

N! = B unc + B nn + B pn + N v (2) 

where N\ is the number of selected events, B unc the un- 
corrected background from natural radioactivity, B nn 



the correlated background from two- neutron captures, 
B pn the correlated background from single-neutron- 
induced events and N u the anti-neutrino signal to be 
measured. We then reverse the selection by imposing 
neutron cuts on the prompt part and positron cuts on 
the delayed part, obtaining: 

N 2 = B' unc + B' nn + e 1 B pn + e 2 N v (3) 

Since both uncorrelated and two-neutron backgrounds 
are symmetric under this selection swap, we have B' unc = 
B unc and B' nn = B nil . Indeed both B' unc and B unc can be 
measured independently and are found to be the same in 
Palo Verde data. The terms B pn and N v are not a priori 
symmetric and we use the factors t\ and e 2 to describe 
the efficiencies for the swapped selection. 

It is essential to realize here that this procedure can 
only be applied if the the trigger system treats the two 
parts of the event in an identical fashion. At Palo Verde 
the symmetric trigger conditions used |]l6f were found to 
have an efficiency very similar to the one that would be 
obtained by separately optimizing the patterns for the 
positron and the neutron parts of the events. 

We can now calculate the difference 

N!-N 2 = (l- ei)B pn + (1 - e 2 )N v (4) 

where B unc and B nn have been eliminated and e 2 can be 
easily obtained from the v c Monte Carlo simulation since 
this process is well known. The derivation of t\ is more 
involved since, as already discussed, the neutron back- 
ground is not easy to model. Here we remark, however, 
that the "swap method" owes its power to the fact that, 
as it will be shown later, t\ ~ 1 and e 2 ~ 0. A small 
(1 — £i) relaxes the accuracy requirements on the Monte 
Carlo simulation to be used to estimate (1 — ei)B pn . 

The Palo Verde detector JIt]] is shown in Fig. All our 
simulations use the Monte Carlo program GEANT H 
to describe the detector and the materials surround- 
ing it. Electromagnetic interactions are simulated by 
GEANT while hadronic interactions are simulated by 
GFLUKA @. Low energy neutron transport is simu- 
lated by GCALOR |po|| . Cuts for tracking neutrons are 
set to 1 MeV for concrete and earth, 100 keV for the 
veto scintillator, 10 keV for the water shielding and 10 -5 
eV for the central detector. Light quenching for protons 
in liquid scintillator is also included pi]. Our program 
successfully simulates the behavior of neutrons from Am- 
Be and e + and 7's from 22 Na sources, which proves that 
neutrino signals are simulated correctly E^j. 

In estimating B pn we consider both the process of 
muon spallation and capture. Each process may happen 
either in the laboratory walls or inside the boundaries 
of the veto system. Our selection cuts for the positron 
part of the v c event have been found to have negligi- 
ble efficiency for muon-induced radioactivity, so that this 
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FIG. 1. Schematic view of the Palo Verde neutrino detec- 
tor. The liquid scintillator is loaded with 0.1% Gd in order to 
reduce the neutron capture time to 30 fis and provide a large 
capture signal (8 MeV 7 cascade). 



phenomenon represents a negligible fraction of the back- 
ground and is not analyzed further. Neutrons from other 
processes such as neutrino interactions with the rock, 
photo-nuclear reactions associated with electromagnetic 
showers generated by muons, and muon elastic scatter- 
ing, are also found to be negligible. 



A. Muon spallation inside laboratory walls 



As it turns out, muon spallation in the concrete walls of 
the laboratory is the dominant component of B pn in the 
case of Palo Verde. Although absolute rate predictions 
are not particularly reliable, in our case a normalization 
point can be obtained from data where the prompt part 
of the event has an energy in excess of 10 MeV. In fact, 
above this energy there is no anti-neutrino signal or neu- 
tron capture but only proton recoils from neutron colli- 
sions. So we use the simulation only to obtain the ratio 



B™°(E < 8 MeV) 
r = B p M n c (B > 10 MeV) 

and then find B pn normalizing to data: 

B pn = r ■ B^ ta (E > 10 MeV). 



(5) 



(6) 



We then determine (1 — e\)r using the Monte Carlo sim- 
ulation. 

As mentioned above, the energy spectrum of neutrons 
from muon spallation is not very well known and a broad 
range of results can be found in the literature. Bar- 
ton [O] suggested that the spectrum of neutrons from 
hadronic cascade follows E~ x / 2 between 10-50 MeV, 
while the spectrum of neutrons from 7r~ capture follows 



a flat spectrum up to 100 MeV. Perkins suggested 
that the neutron spectrum from muon spallation follows 
E~ 16 . The combination of (9.7-E -1 / 2 + 6.0e- B / 10 ) has 
been used in a measurement |j| at a shallow site. It 
has also been suggested Q to use proton and neutron 
spectra following E~~ 1S6 as measured at accelerators for 
photo- nuclear interactions. Finally the Karmen experi- 
ment reported a visible energy spectrum following e ~ B / 39 
for spallation neutrons [jl3| . 

We conservatively choose four spectra, including all the 
options described, as input to our Monte Carlo calcula- 
tion of backgrounds. Tab. | shows the two extreme cases 
of E~ - 5 and E~ 2 - together with the exponential spec- 



trum e 



-E/39 



(Models A). We assume that neutrons are 
produced isotropically. In addition we compute the neu- 
tron spectrum by producing cosmic muons in the energy 
range 0.01 MeV < E^ < 500 GeV according to the proper 
energy and angular distributions f24j . We then gener- 
ate real bremsstrahlung 7's according to the distribution 
l/Ej in the energy range 10 MeV < E 1 < E^. Neutrons 
are then produced from photo-nuclear processes with a 
spectral shape E~ lm Q and an angular distribution 
from H| in the energy range 10 MeV < E n < E 1 . The 
result of this method is also reported in Tab. | as Model 
B. 



Model 



-2.0 



\ A 
A 



E/39 A 
1,86 ' B 



Average 



1.16 ±0.07 
1.20 ±0.11 
1.06 ±0.07 
1.15 ±0.06 



0.69 ± 0.04 
0.67 ±0.07 
0.77 ±0.05 
0.76 ± 0.04 



(l-ei)r 



-0.11 ±0.05 
-0.13 ±0.07 
-0.05 ± 0.06 
-0.11 ±0.04 



-0.10 ±0.05 



BpTtd- 1 ) 
E > 10 MeV 



155 
1.7 
17 
32 



TABLE I. Results of Monte Carlo simulation for neutrons 
produced in the laboratory walls by muon spallation. The 
errors shown are due to limited Monte Carlo statistics. The 
estimated background rate above 10 MeV shown in the last 
column refers to an independent calculation, described in the 
text. The values of B^p(E > 10 MeV) should be compared 
with a total rate of 13.5 ± 0.4 d _1 obtained from Palo Verde 
data. The differences between Models are discussed in the 
text. 



In order to obtain the results in Tab. |[ we use only 
interactions in a 1 m thick concrete shell as neutrons 
produced at larger depths are completely absorbed. The 
10 MeV low-energy cutoff used in the calculations is jus- 
tified by the fact that softer neutrons are completely ab- 
sorbed by the 1 m thick water buffer surrounding the 
Palo Verde central detector. From Tab. | we can see 
that the proton-recoil energy is only weakly dependent 
upon the neutron energy so that both e\ and r remain 
almost constant for drastically different neutron energy 
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spectra. Furthermore, the uncertainties on e\ and r have 
little effect on the factor (1 — e\)r that directly enters the 
neutrino measurement. This implies that the neutron- 
capture signal, common to both the neutrino signal and 
the B pn background, is similar to the proton-recoil signal 
of the background, but different from the positron signal 
of a neutrino event. 

In Fig. 2 we show the energy deposited by the neutron- 
induced proton recoil in the most energetic cell of the 
prompt part of the events. The four different neutron 
spectra used are normalized to data for energies above 
10 MeV. 

To verify the results obtained we independently calcu- 
late the spallation background by using a neutron yield 
of 6 x /i g cm for normalization. This number 

is obtained by rescaling the measurements of ^ to our 
depth of 32 m.w.e. A total of 7.8 x 10 6 neutrons are gen- 
erated daily in our lab walls. The calculated background 
rates B^(E > 10 MeV) are shown in the last column 
of Tab. |, which can be compared with our measurement 
B^ ta {E > 10 MeV) = 13.5 ± 0A(stat.) d" 1 . It is clear 
that our measurement falls somewhere in the middle of 
the predictions and the spectra chosen for the simulation 
cover a conservative range of possibilities. We conser- 
vatively maintain all four options and use the possible 
differences as contributions to the systematic errors. 

Finally we average (1 — e\)r from Tab. [j| obtaining 
—0.10 ± 0.05, and then proceed to calculate the back- 
ground from neutron spallation in the walls as (1 — 
ei )B pn = 1.35 ±0.68 d _1 . 



Model 



£-"- b , A 
E- 2 - , A 
e"^ 39 , A 

E -1M B 



El 



1.17 ±0.12 
0.94 ±0.08 
0.97 ± 0.04 
1.13 ±0.05 



0.72 ± 0.07 
1.41 ±0.15 
1.12 ±0.05 
1.12 ±0.05 



(1 



ei ir 



-0.12 ±0.08 
0.08 ±0.11 
0.03 ± 0.04 

-0.15 ±0.06 



B, 



(d- 1 ) 
E > 10 MeV 



2.2 
0.06 
0.8 
0.9 



TABLE II. Results of Monte Carlo simulation for neutrons 
produced in the water buffer by muon spallation. The errors 
shown are due to limited Monte Carlo statistics. The esti- 
mated rate above 10 MeV shown in the last column refers 
to an independent calculation, described in the text. It is 
clear that the rates found for this channel are negligible with 
respect to the rates in Tab. Ill 



(compared to spallation), with an upper limit around 
100 MeV (muon mass). The underground laboratory 
at Palo Verde is built with low activity concrete using 
dolomite as aggregate. The elemental composition of 



concrete is shown in Tab. HI together with the muon cap- 
ture rate and the neutron yield per capture for each ele- 
ment. Almost every capture produces one neutron. The 
total un- vetoed muon rate in the walls is 2 kHz, of which 
0.9 kHz is due to fx~ . The stopping fjf rate amounts to 
90 Hz. Using Tab. Ill we obtain a total muon capture 
probability of 67%, resulting a neutron production rate 
of 60 Hz in the laboratory walls. 



B. Muon spallation inside the veto system 

Inefficiencies of the cosmic-ray veto system result in a 
component of B pn from neutrons produced within the de- 
tector. The veto inefficiency at Palo Verde is measured to 
be (0.07 ± 0.02)% for through-going muons (two missed 
hits). Only neutrons produced in the water buffer are 
important here since muons responsible for neutron spal- 
lation directly in the central detector scintillator would 
be easily detected and discarded. 

Using the same procedure as above we obtain a neutron 
yield of 1600 d _1 from the water buffer. The correspond- 
ing ei, r and (1 — ei)r are given the in Tab. || together 
with background estimates B^°(E > 10 MeV). It can be 
readily seen from the Table that B^ C (E > 10 MeV) in 
the water buffer has a negligible rate compared to that in 
Tab. |, and their (1 — ei)r are very similar. Hence in the 
rest of our calculations we will neglect this contribution. 



C. Muon capture inside laboratory walls 

The muon capture process is rather well understood 
and the resulting neutrons tend to have a soft spectrum 



Element 



H 

C 
O 
Al 
Si 
Mg 
Ca 
Fe 



Fraction by 
mass (%) 



0.6 
10.4 
50.6 
0.3 
1.2 
10.7 
22.9 
3.3 



Capture rate 
(10 5 s" 1 ) 



0.004 
0.388 
1.026 
7.054 
8.712 
10.67 
25.57 
44.11 



[ 25| 
[ 25| 
[25| 
[ 25| 
[ 25| 
[ 25| 
[25| 
[25 1 



n yield/capture 
(/i-capture -1 ) 



26 



1 

~ 1 
0.98 
1.26 J27| 
0.86 |27 

1 * 
0.75 §7| 
1.12 [27] 



*actual value not known, assumed to be 1 

TABLE III. Elemental composition together with muon 
capture rates and neutron yields for the Palo Verde concrete. 
The concrete contains 3% reinforcing steel, 16% cement and 
81% crushed dolomite aggregate. 



The neutron energy spectrum from capture is simu- 
lated taking into account both the soft neutrons from 
nuclear evaporation and the hard neutrons from direct 
emission. For light elements such as 12 C and 16 0, indi- 
vidual lines are present in the neutron spectrum p8| , p9| , 
while for heavier elements such spectrum has the proper- 
ties of a continuum. We use the energy spectra in ref. |28] 
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FIG. 2. Energy deposited in the most energetic cell for different neutron spectra. Data and Monte Carlo total rates are 
normalized above 10 MeV where there is no i> e signal in the data. The dashed line is the MC prediction for E < 10 MeV. 



to simulate neutrons from 12 C and 16 O, while those 
from heavy elements are simulated according to [p7| . 
Capture on hydrogen happens at a negligible rate and 
is disregarded here. From the simulation, we obtain 
B pn = 0.10 ± 0.05 d" 1 and (1 - e x ) = 0.23 ± 0.32 so that 
(1 - ei)B pn = 0.02 ± 0.03 d" 1 where the error includes 
Monte Carlo statistics and all systematic uncertainties. 
We conclude that this background is negligible compared 
to other channels. 

D. Muon capture inside the veto system 

In analogy to the spallation case, neutrons from muon 
capture inside the veto system can contribute to B pn for 
un-tagged muons. Even considering the conclusions from 



the previous two sections this background cannot be a- 
priori dismissed as negligible since the veto counter inef- 
ficiency for single hits (such as would result from a stop- 
ping muon) is measured to be (4 ± 1)%. In the buffer- 
water muon capture on 16 O is the only significant source 
of neutrons since the capture rate on hydrogen is very 
small. 

The total fi~ rate in our detector is about 860 Hz, 
of which 86 Hz are stopping muons. This results in a 
rate of un-tagged neutrons of 52500 d _1 . Using the en- 
ergy spectrum from ref . ^8| , we obtain from Monte Carlo 
simulation B pn = S^iO^d" 1 and (1-ei) = 0.22±0.03. 
Finally, this background contributes (1 — ei)B pn = 0.86± 
0.50 d _1 where, as usual, the error includes all uncertain- 
ties. 
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Rate (d _1 ) 


(1 — ei)-Bpn (Spallation in walls) 
(1 — £i)-Bpn (Spallation in detector) 
(1 — ei)-B pll (Capture in walls) 
(1 — £i)^?pn (Capture in detector) 


0.19 ±0.26 
— 08 ± 08 


Total (1 -ei)B P n (MC) 


0.11 ± 0.27 


Ni (Data) 
N 2 (Data) 
Ni - N 2 (Data) 


8.75 ± 0.28 
9.07 ±0.29 
-0.32 ±0.20 



TABLE IV. Comparison of data and Monte Carlo for an 
event selection with no efficiency for the anti-neutrino signal 
(see text). Errors are statistical only. Note that Ni and N 2 
are correlated. 



E. Verification of the method 

In summary, all the above backgrounds add to a total 
rate (1 - ex)B pn = 0.5 ± 0.8 d -1 , very close to 0. The 
error is dominated by systematics, particularly stemming 
from uncertainties in the neutron energy spectrum. 

In order to verify the correctness of the method, we 
can directly measure in the data a similar background 
by slightly modifying the anti-neutrino selection cuts so 
that no signal is detected. Positrons (from z7 e interac- 
tions) differ from proton recoils (from background neu- 
trons interactions) by the annihilation 7's with energies 
of less than 511 keV. An event selection requiring more 
than 600 keV for each of the hits will result in the total 
rejection of the neutrino signal. Hence in this case only 
the -Bp n term will be present after swap selection: 



JVi - N 2 = (1 - ei) • B P 



(7) 



Tab. 
only. 



IV 



shows the result of this test with background 
The values of N\ — N2 from data is consistent, 
within errors, with Monte Carlo estimate of (1 — ei)-B pn . 
Different selection cuts for the background yield similar 
results. 



IV. COMPARISON WITH "ON-OFF" 
BACKGROUND SUBTRACTION 

The advantages of the swap method become clear in 
the comparison with the "ON-OFF" method, as shown 
111 Tab. [V| The Table summarizes the results of the 
Palo Verde experiment for the 1999 data taking period 
from |Q. Ei is indeed very close to 1 resulting in a very 
small residual background (1 — ei)B pn , even for a back- 
ground rate B UDC + B nn + B pn as high as 27 d _1 . In the 
new variable iVi — N%, not only the terms B unc and B nn 
drop completely, but also B pn is strongly suppressed. A 
conservative 160% uncertainty on (1 — ei)B pn only cor- 
responds to a 4% error on N„. On the other hand, €2 is 





1999 "ON" 


1999 "OFF" 


No. of days 
i/f, efficiency 

£2 


110.95 
0.112 
0.159 


23.40 
0.111 
0.159 


(1 — £l)-Bpn, spall, in walls (d i ) 
(1 £l)^?pn, spall, in water (d ) 
(1 — £l)5 pnj C apt. in walls (d 1 ) 
(1 £l)-^pn, capt. in water (d ) 


-1.35 ±0.68 

0.02 ± 0.03 
0.86 ± 0.43 


-1.33 ±0.67 

0.02 ±0.03 
0.86 ± 0.43 


Ni (d" 1 ) 
No (d' 1 ) 
N„ (d" 1 ) 


52.9 ±0.7 
32 3 ± 5 
25.2 ± 0.9 


43.9 ± 1.4 
31.7 ± 1.2 
15.1 ± 1.9 


Bunc + Bnn + Bpn (d X ) 


27.7 ±0.6 


28.8 ± 1.3 


u e observed (d _1 ) 
v e expected (d _1 ) 


225 ±8 
218 


136 ±17 
130 



TABLE V. Palo Verde results from 1999 data taking. Er- 
rors are statistical except for (1 — ei)_B pn where errors are sys- 
tematic. The individual background rates are approximately 
4 d _1 for B U nc, 14 d" 1 for B nn and 10 d" 1 for B pn 



only 0.16, small enough that the statistical power of the 
D e signal is retained. 

Correcting Ni in both columns by their respective ef- 
ficiencies and subtracting column 2 from column 1 in 
the Table, we find that the "ON-OFF" method gives a 
neutrino rate of 77 ± 14(stat.) ± 8(syst.) d _1 for an ex- 
pectation of 88 d _1 in the no-oscillation hypothesis. In 
calculating the signal essentially only one reactor out of 
three is used, while the statistical fluctuations in the flux 
of all reactors along with the background contribute to 
the errors. The systematic error includes uncertainties 
on positron and neutron efficiencies (5%), v e selection 
(8%) and v e flux estimate (3%). 

In the case of the swap method, we find 225±8(stat.)± 
17(sys.) d" 1 (137±17(stat.)±14(sys.) d" 1 ) for high (low) 
power against a prediction of 218 d _1 (130 d _1 ) for the 
case of no oscillations. Here all reactors contribute to the 
signal and, in fact, the contributions for the two periods 
with different power can be used together to strengthen 
the measurement. Indeed the statistical error drops from 
18% in the case of "ON-OFF" to 3.5% (12%) for high 
(low) power. While systematic errors from efficiencies 
and flux are the same as in the previous case, the error on 
D e selection is now only 4% because some of the selection 
systematics cancel in the N\ — N2 difference. A new 
uncertainty due to the B pn estimate (4%) appears. 



V. SUMMARY AND DISCUSSIONS 

We have shown that a novel method of analysis, ap- 
plicable to low energy neutrino experiments using cor- 
related signatures, provides substantially more accurate 
background subtraction over more traditional techniques. 
While the new method was applied first to a reactor 
neutrino oscillation experiment, it can be more generally 
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used in experiments where: a) the signal events consist 
of two sub-events correlated in time or space, b) the two 
sub-events are distinctively different from each other in 
signal but similar in backgrounds or vice- versa, and c) the 
detector and trigger treat the two sub-events in identical 
fashion. 

These criteria apply to several future neutrino experi- 
ments such as KamLAND and LENS. In KamLAND |] 
electron anti-neutrinos from reactors will be detected in 
1 kton liquid scintillator as a positron with energy de- 
posit between 1 and 8 MeV correlated in time with a 
neutron which gives a 2.2 MeV 7 line from capture on 
protons. The correlated neutron background, which in- 
cludes both one- and two-neutron events, is expected to 
have the same magnitude as the random background. 
The application of the method is therefore straightfor- 
ward: both random and two-neutron backgrounds can 
be eliminated and the one-neutron background can be 
estimated in a way which is very similar to what we 
discussed above. The LENS experiment || is designed 
to detect solar neutrinos via, for example, the process 



160 



Gd 



160 



Tb* , where the signature consists of 



an electron (0.04-2 MeV) and a 7 (64 keV) correlated in 
time. One of the main backgrounds is the random coin- 
cidence of 7's from radioactive impurities in the detector 
and it can be easily suppressed by the method described 
above. 
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